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Abstract —In this paper, we consider a down-link transmission 
of multicell virtualized wireless networks (VWNs) where users 
of different service providers (slices) within a specific region 
are served by a set of base stations (BSs) through orthogonal 
frequency division multiple access (OFDMA). In particular, we 
develop a joint BS assignment, sub-carrier and power allocation 
algorithm to maximize the network throughput, while satisfying 
the minimum required rate of each slice. Under the assumption 
that each user at each transmission instance can connect to no 
more than one BS, we introduce the user-association factor (UAF) 
to represent the joint sub-carrier and BS assignment as the opti¬ 
mization variable vector in the mathematical problem formula¬ 
tion. Sub-carrier reuse is allowed in different cells, but not within 
one cell. As the proposed optimization problem is inherently 
non-convex and NP-hard, by applying the successive convex ap¬ 
proximation (SCA) and complementary geometric programming 
(CGP), we develop an efficient two-step iterative approach with 
low computational complexity to solve the proposed problem. 
For a given power-allocation, Step 1 derives the optimum user- 
association and subsequently, for an obtained user-association, 
Step 2 find the optimum power-allocation. Simulation results 
demonstrate that the proposed Iterative algorithm outperforms 
the traditional approach in which each user is assigned to the BS 
with the largest average value of signal strength, and then, joint 
sub-carrier and power allocation is obtained for the assigned 
users of each cell. Especially, for the cell-edge users, simulation 
results reveal a coverage Improvement up to 57 % and 71 % for 
uniform and non-uniform users distribution, respectively leading 
to more reliable transmission and higher spectrum efficiency for 
VWN. 

Index Terms —Complementary geometric programming, suc¬ 
cessive convex approximation, joint user association and resource 
allocation, virtualized wireless networks. 

I. Introduction 

A. Motivation 

The next generation of wireless networks encounter dif¬ 
ferent types of challenges such as meeting the requirements 
of new high-data rate applications, reducing the cost, and 
increasing the spectrum and infrastructure efficiency. In this 
context, recently, virtualization of wireless networks has been 
brought up as a promising approach to tackle these issues. In 
virtualized wireless networks (VWNs), the physical resources 
of one network provider, e.g., spectrum, power, and infrastruc¬ 
ture, are shared between different service providers, also called 
slices m, m, 0. Generally, each slice comprises of a set of 
users, and has its own quality-of-service (QoS) requirements. 
To harvest the potential advantages of VWN, effective resource 


allocation is a major concern which has been drawing a lot of 
attention in recent years. 

For instance, in m, a resource management scheme is 
studied by introducing two types of slices, including rate- 
based and resource-based slices, where the minimum rate 
and minimum network resources are preserved for each slice, 
respectively. Furthermore, in a, interactions among slices, 
network operator, and users are modeled as an auction game 
where the network operator is responsible for spectrum man¬ 
agement on a higher level, and each slice focuses on QoS 
management for its own users. To preserve the QoS of slices 
from wireless channel fading, the admission control policy 
is proposed in a where the requirement of each slice is 
adjusted by the channel state information (CSI) of its own 
users. To extend the feasibility condition of VWN in order 
to support diverse QoS requirements, 0 considers the use 
of massive MIMO where the access point is equipped with 
a large number of antennas. In Q, the combination of time, 
space and elastic resource allocation for OFDMA system is 
considered. Advantages of full-duplex transmission relay in 
VWN have been investigated in jS]. 

Generally, these works have focused on analyzing the 
resource allocation problem in a single-cell VWN scenario. 
However, in practice, the coverage of a specific region may 
require a set of BS, in a multi-cell VWN scenario. The key 
question in a multi-cell VWN scenario is how to allocate the 
resources to maintain the QoS of each slice, while improving 
the total performance of VWN over a specific region. In this 
paper, we consider a multi-cell OFDM based VWN where 
the coverage of a specific region is provided by a set of BSs 
serving different groups of users belonging to different slices. 
The QoS of each slice is represented by its minimum reserved 
rate. Each user of each slice can be only served by one BS and 
this BS is not predetermined by the distance or by measuring 
the average received signal strength of BSs. Consequently, in 
this setup, the resource sets in the related optimization problem 
involve the set of BSs, sub-carriers and power for each user 
belonging to each slice. 

In this paper, the objective of proposed resource allocation 
problem is to maximize the total throughput of VWN in the 
specific region subject to power limitation of BSs, minimum 
required rate of each slice, and sub-carrier and BS assignment 
limitations. Based on the limitations of down-link OFDMA 
transmission, each sub-carrier can be assigned to one user 
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within a cell and each user can be associated to only one 
BS. Since in this optimization problem, the sub-carrier as¬ 
signment and BS association are inter-related, we introduce 
the user association factor (UAF) that jointly determines the 
BS assignment and sub-carrier allocation as the optimization 
variable vector. Due to this user-association constraint and the 
inter-cell interference, the proposed optimization problem is 
non-convex and NP-hard, suffering from high computational 
complexity E). We apply the frameworks of complementary 
geometric programming (CGP) and the successive convex ap¬ 
proximation (SCA) approach ifTOl . ifTTl . ifT^ . ifTSl to develop 
an efficient, iterative, two-step algorithm to solve the proposed 
problem. For a given power-allocation. Step 1 derives the 
optimum user-association solution, and subsequently, with this 
obtained user-association solution. Step 2 finds the optimal 
power allocation. Such 2-step optimization process is repeated 
until convergence. Furthermore, it can be shown that even 
the simplified problem of each step still involves non-convex 
optimization problem. Thus, by applying various transforma¬ 
tion and convexification techniques, we develop the analytical 
framework to transform the non-convex optimization problems 
encountered in each step into the equivalent lower-bound 
geometric programming (GP) problems, ca, which can be 
solved by efficient softwares, e.g., CVX IITtII . 

Simulation results demonstrate that the proposed approach 
can significantly outperform the traditional scenario where 
the BS assignment is based on the largest average signal- 
to-interference-plus-noise ratio (SINR), and subsequent sub¬ 
carrier and power allocation is derived for the users of each 
cell. The simulation results reveal that considering UAF can 
increase the feasibility of resource allocation problem (the 
required rate of each slice will be satisfied with a higher 
probability as compared to the traditional approach). Specif¬ 
ically, the proposed algorithm can significantly increase the 
probability of achieving higher rates for the cell-edge users, 
resulting better coverage for the VWN. 

B. Related Works 

Our work in this paper lies along the intersection of two 
research contexts in resource allocation problems: 1) multi¬ 
cell OFDMA wireless networks, and 2) VWNs. 

There exists a large body of research conducted in re¬ 
source allocation for multi-cell OFDMA wireless networks. 
For example, in IITSl . the resource allocation in conventional 
OFDMA-based network is studied using BS-assignment based 
on the largest average received signal strength from the BS at 
each user. An iterative algorithm for maximizing the weighted 
sum of minimal user rates in each BS is explored in ifTTll . Joint 
cell, channel and power allocation in multi-cell relay networks 
is explored in ifTSl . where each user is assigned to the BS with 
the highest channel gain. In 1191 . a proportional fair resource 
allocation in a multi-cell OFDMA networks is proposed with 
the aim of maintaining the quality of experience of users by 

'This average is derived based on the measurement of users over one 
specific window in both idle and active phases where the size of window and 
the measurement of each user are adjusted based on specification of wireless 
network standards na. 


considering a utility function based on the mean opinion score. 
In II 20 I . joint scheduling of resource blocks, power alloca¬ 
tion, and modulation and coding scheme in LTE-A system 
is considered by using the criteria of proportional fairness. 
A similar problem in OFDMA cognitive radio networks is 
studied in m] where an iterative algorithm is proposed to 
solve the sub-carrier and power allocation. Similarly, in ll22ll . 
a resource allocation problem for jointly optimizing the energy 
and spectral efficiency is proposed for a multi-cell OFDMA 
wireless network by considering an energy and spectral effi¬ 
ciency trade-off metric. In ||23]| . the authors have considered 
an energy efficient resource allocation problem for a multi-cell 
OFDMA network in a conventional wireless network where 
the available values of channel state information (CSI) are 
imperfect. In ll24ll . a resource allocation algorithm is proposed 
for a two-cell down-link OFDMA network with a fractional 
frequency reuse scheme between BSs. 

In the aforementioned works (i.e., Qa, ini, m, ma, 
Eol, ED, Ea, (ED), the BS assignment algorithm is sepa¬ 
rated from the sub-carrier allocation, while joint sub-carrier 
and power allocation is applied for multi-cell scenario. Com¬ 
pared to this research direction, we consider UAF which 
jointly assigns the BS and sub-carrier for each user. In our 
proposed algorithm, after assigning UAF, the power allocation 
algorithm is applied. Furthermore, we maintain the implemen¬ 
tation limitations of multi-cell OFDMA networks by proposing 
new constraints. 

As previously mentioned, the resource allocation in VWNs 
has received growing attention. In different aspects of 
VWN including resource discovery and allocation as well 
as the research challenges have been discussed. Besides m, 
a, 0 , a, Q, in E5I . the challenge in allocating physical 
resource blocks (PRBs) to various slices in an LTE network 
has been addressed considering a single BS scenario. In EH, 
an opportunistic algorithm to allocate the resources to virtual 
operators is proposed by differentiating the resource require¬ 
ments among operators as baseline and fluctuate requirements 
to ensures the minimum QoS requirements of each virtual 
operator. In El^ the concept of virtualization has been extended 
to a LTE network by considering virtual operators or slices 
each with various bandwidth requirements in terms of the 
physical resource blocks (PRBs) in LTE. To the best of our 
knowledge, the multi-cell scenario of VWNs has not been 
studied in the previous related works in this context. 

C. Structure of Paper 

The rest of this paper is organized as follows. In Section 
II, the system model and problem formulations are introduced. 
Section III contains mathematical background, the detailed ex¬ 
planation of iterative algorithms and computational complexity 
analysis. Section IV demonstrates the simulation results and 
their detailed analysis followed by concluding remarks in 
Section V. 

II. System Model 

We consider the down-link transmission of a VWN where 
the coverage of a specific area is provided by the set of BSs, 
i.e., A4 = {1,..., M}. The total bandwidth of B Hz is divided 
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Fig. 1; Illustration of multi-cell VWN with M BSs serving 
sets of users belonged to different slices. 


into a set of sub-carriers, 1C = ^K} and shared by 

all BSs through orthogonal frequency-division multiple-access 
(OFDMA). The bandwidth of each sub-carrier, i.e., 
is assumed to be much less than the coherent bandwidth of the 
wireless channel, meaning that the channel response in each 
sub-carrier is flat. This set of BSs serves a set of slices, i.e., 
Q = {1,..., G}, where the slice g has a set of users (denoted 
by Afg = {1,..., Ng}) and requests for a minimum reserved 
rate of A typical illustration of this setup is shown in 
Fig. [U where N = Thg^g *^he total number of users. 

Let hm,k,ng and Pm,k,ng be the channel state information 
(CSI) of the link from BS m G Af to user Ug of slice g 
on sub-carrier k and allocated power of user Ug of slice g on 
sub-carrier k, respectively. In this scenario, due to the OFDMA 
limitation, each user is assigned to one BS, and to avoid intra¬ 
cell interference, orthogonal sub-carrier assignment is assumed 
among users in a cell. We assume that there is no pre-allocated 
BS for users, and consider /3m,fe,ng G {0,1}, which represents 
both sub-carrier allocation and BS assignment indicator for 
user Ug of slice g on sub-carrier A: of BS m which has been 
called user association factor (UAF). Note that I3m,k,n = 1 
when BS m allocates sub-carrier k to user Ug, and I3m,k,n = 
0, otherwise. 

Consider P = ^ = 

[Pm,k,ng]y^g^^ as the vector of all transmit powers 
and UAFs of users, respectively. Therefore, the rate of user 
Ug at sub-carrier fc of BS m can be expressed as 


Rm,k,ng (P) — lo§2 


1 + 


Pm,k,rig hm,k,ng 


( 1 ) 


where Im,k.ng = 

g^g^^Vn^Gj^fg ,n^g^ng Rm' ,k,n'g^m,k,ng tS 

the interference observed by user Ug in cell m and sub-carrier 
k, and is the noise power. Without loss of generality, noise 
power is assumed to be equal for all users in all sub-carriers 
and BSs. Now, from ([T]i, the required minimum rate of slice 
g & G can be represented as 

Cl • E EE P^PngRm,k,ng{V) > R^ ^9 G G- 

mUM rigGAfg k&K 

Considering the maximum transmit power limitation of each 
BS, 

C2 : E E E Prn,k, Ug < PTi Vm G M, 

gsp rig&N'g k&K 

where is the maximum transmit power of BS m. Fur¬ 
thermore, the OFDMA exclusive sub-carrier allocation within 
each cell m can be expressed as 

C3 : E E Pm, k,Tig <1, Vm G M, Vfc G 1C. 

g&Q rig&Mg 

In this setup, due to the limitation of multi-cell OFDMA, we 
restrict the access of each user by the following constraint 

C4 : [ ^ Pm,k,ng\ [ E E P'T^'A,ng\ = 0, 

k^K ym'^m k^K. 

'iug € Afg, Vg G G, Vm G A4. 

C4 imposes that when sub-carrier k from BS m is assigned 
to user Ug, it will not be assigned to user Ug by BS m' 
where m' ^ m and Mm',m G Ad. C4 is novel since it 
jointly assigns the BS and sub-carrier over the specific region. 
Therefore, we call Pm,k,n as the user association factor (UAF) 
which simultaneously determines the BS and sub-carrier in the 
VWN. The joint power, sub-carrier and BS assignment can be 
formulated as 

^ ^ ^ ^ ^ ^ Pm,k,ngRm,k,ng{P), (2) 

m^M gi^g rigGAfg k(^K 

subject to; Cl - C4. 

The optimization problem introduced in (|2]l has a non-convex 
objective function due to inter-cell interference and involves 
non-linear constraints with combination of continuous and 
binary variables, i.e., P and /3. In other words, © is a non- 
convex mixed-integer, NP-hard optimization problem Et), 0. 
Therefore, proposing the efficient algorithm with reasonable 
computational complexity is desirable. 


III. Two-Level Iterative Algorithm for Joint User 
Association and Resource Allocation 

To tackle the computational complexity of ©, we adopt 
an iterative approach to iteratively And the UAF and power 
allocation for each user in two steps as shown in Algorithm 
1. In Step 1, for given power allocation vector, the UAF 
considered as the variable of the user-association problem is 
solved by Algorithm LA to be discussed in detail in Section 
III.B. 
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Algorithm 1 Iterative Joint User Association and Power 
Allocation Algorithm 

Initialization: Set t := 1, (3{t = 1) = [1], where 1 is a 
vector, and the power over of each sub-carrier of BS 

m is equal to /K. 

Repeat 

Step 1: Derive (3*{t) from (fTsT i by using Algorithm l.A; 

Step 2: For derived f3*(t), find p*(f) from (l2?t by 
applying Algorithm l.B; 

Step 3: Stop if ||/3*(f) — /3*(f — 1)|| < Si, and ||p*(f) — 
p*(< — 1)11 < £ 2 , otherwise, set f := f -f 1 and go to Step 1. 


This derived UAF is then used in Step 2 to find the 
corresponding allocated power as the solution of the power- 
allocation optimization problem by Algorithm l.B to be dis¬ 
cussed in detail in Section II.C. Steps 1 and 2 are iteratively 
executed until both the current UAF and power allocation 
vector solutions are not much different from their values 
obtained in the previous iteration. In other words, the sequence 
of the UAF and power allocation vector solutions can be 
expressed as 

/3(0) ^ P(0) ^ ... I3{t)* P{t)* (3* P* , (3) 

Initialization Iteration t Optimal solution 

where f > 0 is the iteration number and /3(f)* and P(f)* are 
the optimal values at the iteration t from convex transformation 
of related optimization problems in each step. The iterative 
procedure is stopped when 

||/3*(f)-/3*(f-l)|| <£i and ||p*(i) - p*(t - 1)|| < 

where 0 < ei <C 1 and 0 < £2 1. 

Notably, both the user-association and power-allocation 
optimization problems are still non-convex and suffer high 
computational complexity. To solve them efficiently, we apply 
complementary geometric programming (CGP) for each step 
ini in which via different transformations and convexification 
approaches, the sequence of lower bound GP based approxi¬ 
mation of relative optimization problem is solved as described 
in detail in the following sections. 


A. Complementary geometric programming (CGP):A brief 
review 


The standard form of GP is defined as 

min/o(x), (4) 

X 

subject to; 

/.(x)<l, Vz = 0,l,...,/, 

Pj(x) = l, Vj = 0,1,..., J, 

where x = [xi,X 2 , ■■■,Xn] is a non-negative optimization 
variable vector, gj{x) for all j is a monomial function, i.e., 

N 

5j(x) = n CjnXn"' where > 0, ajn e 3?, and /o(x) 

n—1 

and /i(x) for all i are posynomial functions, i.e., /i(x) = 
Sfc=i n Ciknxf’’"- ■ However, in (HJi, there are a lot of restric- 

n—1 

tions on the equality and inequality constraints which cannot 
be met for many practical problems related to the resource 
allocation of wireless networks such as the optimization prob¬ 
lem considered in this paper. For example, in some cases, the 
equality constraints contain posynomial functions, inequality 
constraints present lower bound of posynomial function or the 
posynomial functions contain negative coefficients. Depending 
on the nature of the optimization problem, these types of 
problems belong to either one of classes of optimization 
problems such as generalized GP, signomial programming or 
complementary geometric programming (CGP). A CGP can 
be presented as 

minFo(x), (5) 

X 

subject to: 

F, (x)<l, Vi = l,...,J, 

G, (x) = l, Vj = l,...,J, 

where Uo(x) = f+{x) - ff (x), F)(x) = for all i = 

I,-- - ,I and Gj{x) = in which /o^(x), ff(x) for all 
/ = 0,1, • • • , are posynomial functions, and gj{x) and fj(x) 
are monomial and posynomial functions for all j = 1, • • • J, 
respectively, ||3^ . 

Toward solving 0, one approach is to convert the CGP 
optimization problem into a sequence of standard GP prob¬ 
lems na that can be solved to reach a global solution. 
In other words, successive convex approximation (SCA) is 
applied ll^ . Il34l . where the non-convex optimization prob¬ 
lem is approximated as a convex problem in each iteration. 
Lets consider I as an iteration number for this sequence. 
More specifically, arithmetic-geometric mean approximation 
(AGMA) can be applied to transform the non-posynomial 
functions to posynomial form, i.e., Fi(x), and Gj{x) to its 
monomial functions, respectively. 


Geometric programming (GP) is a class of non-linear 
optimization problems which can be solved very efficiently 
via numerical methods HD. Thus, a significant amount of 
research has been done in order to convert the resource 
allocation problems into GP problems, so that it becomes 
computationally tractable ifTOll . ifTTll . 1^ . 1^ . ll3D . 


Using AGMA, at the iteration I, the approximated forms of 
/r(x) = Ef=i 9l (x) and /^ (x) = dli^) are 


fi (x(o)= n 


K.- 

TT 1 


11 1 

k^l 

[ <~ii) J 


“fc (0 


(6) 
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and 


/j(x(0) 


TT 1 

(gim)\ 

11 1 

K Ciii) ) 


can 


(7) 


where {1) 


9l (’‘(i-1)) 


and Cfc(0 


ffiwi-i)) 


. Now, 


we have Fi(x{l)) = ^ — and G,-(x(t)) = 

which are posynomial and monomial functions, respectively 
112. Now, the optimization problem related to each iteration 
Z of (|5]l is 


niin S + f+ (x(0) - /o (x(0), (8) 

x(t) 

subject to: 

G,(x(0) = l, Vj = l,...,J, 

where S > 0 is a sufficiently large constant which is added to 
the objective optimization problem ([8]l to ensure that the ob¬ 
jective function is always positive El. However, the objective 
function of (HJ still cannot satisfy the posynomial condition of 
(Hi. To take the hnal step to reach the GP based formulation for 
each iteration, we introduce the auxiliary variable a:o > 0 to 
reach a linear objective function and apply it for transformation 
of ® into 

mina:o(0, (9) 

Xo(t) 

subject to: 

S + /o+(x(0) ^ 

/o" (x(0) + a:o “ 

F, (x(f))<l, Vi = 0,1,...,/, 

G, (x(f)) = l, Vj = l,...,J, 

where Xo(f) = [xo{l),Xn{l), - ■ • ,xo(0]. Similar to /j(x), 

~+ (xfZ)) 

term can be converted into posynmial function via 

/o WO)+ 3:0 

AGMA, and hnally, the resulting optimization problem has a 
GP based structure and can be solved by efficient numerical 
algorithm or available softwares lEl- 

It has been shown that the optimal solution of iterative 
algorithm based on the approximation of problem (|5]l into its 
GP based approximation has a very close performance to the 
optimal solution ITT^ . 


B. User Association Problem 

Considering a hxed P(f), we have the following optimiza¬ 
tion problem, referred to as the user association optimization 
problem, 

max EE E E /3m.fe.ng/?m.fe.ng(P(f)), (10) 

^ maM gdQ UgGAfg kalC 

subject to: C1,C3, C4, 


where Rm.k,n (P(i)) a hxed value derived based on P(f) 
and 

Cl : E E E Pm,k,ngR'm,k,ng{^{i'^^ ^ Rg 5 

mCiM ngdAfg kalC 

yg&G- 


In ([Tol l, the only optimization variable is (3, and therefore, 
(doll has less computational complexity compared to 0. 
However, it still suffers from the integer optimization variable 
(3. In addition, due to C4 and the objective function, (fT^ 
is still a non-convex optimization problem. To overcome the 
computational complexity of (fTOl i. we hrst relax the UAF 
variable as ^rn,k,n S [O;!]- Also, we apply the technique 
in Section III. A to convert (fTol l into the GP formulation. 

To reach to a standard GP formulation, the equality con¬ 
straint in C4 should involve only monomial functions. In the 
following, we hrst relax and then apply iterative AGMA algo¬ 
rithm (as in 0 and 0) to get the monomial approximation for 
C4. Also, we show how we can convert the objective function 
of 0 into the standard form of GP. 

Proposition 1: Consider ti as the index of iteration for 
solving ([Toll, Xjn,ng{tl) = Y.k&K.^rn,k,ng{.ti) and = 

Y.meMT.k&ic l^rn,k,ng{ti). C4 Can be approximated by the 
following constraints. 


C4.1: 


C4.2: 


C4.3: 


C4.4: 


^m,ngi^t) 3- Xm,ng(ti)yng(ti)Sg^ ,i^^{ti) < 1 

Vrig e Afg, g € G, Vm G M, 

r 1 


t^m,ng (G) ^ 


,(*i 


< 1 , 


Ani^rig (G ) 

xl^,ngitl) 

_ Oim,ng (^1) _ 

'irig G Afg, g G Q, Vm G M, 

(, \ TT (G)"' ^ 

,(G) [[ 


kGK 


^m,k,ng (^1) 

Vn^ € Afg, g gQ, Vm G A4, 


Pm,k,rig (^1) 


= 1 , 


yngih) n 

mGAi,k^lC 


n -Vni,k,ng(tl) 


'kjm.k.ng (G) 


= 1 , 


Vrig G Afg, g € G, Vm G A4, 
where Sm,rtg(G) is an auxiliary variable, and. 


Am,rig (G) — 

0^m,ng (G ) — 


1 


k'm,k,ng (G ) — 

gm,k,ng (G) 


Xm,ng{h - 1) + 1’ 
Xm,ng (^1 i^) 

a;m,ng(G - 1) + 1’ 

i3m,k,ng{il 1) 


J2keK l^^,k,ng{ti 1) 

Pm,k,ng{Al 1) 

J2mGM '^k^K Prn,k,ng (G — 1) 
Vrig G Afg, g &G, yrn G A4. 


( 11 ) 

( 12 ) 

(13) 

(14) 


Proof: See Appendix A. ■ 

Based on C4.1-C4.4, C4 is transformed and represented by 
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Algorithm l.A ; Algorithm to derive User Association Factor 
Initialization: Set ti := 1, (3{ti = 1) = j3{t) and P(<i = 
1) = P(t:) and set arbitrary initial for Sm,ng{ti), 

Repeat 

Step 1- Update ' 7 m,ng(U )5 Q^m,ng(U); 

Pm,fc,ng(fi)i Cm,fc,ng(U)? and from (I1 llt- 

(0 and ([I7]l-([l9ll, 

Step 2: Find optimal UAF in (flST l using CVX, 


a locally optimal solution that satisfies the KKT conditions of 
the original problem. 

Proof: In HI, it is shown that the conditions for the 
convergence of the SCA lf3^ are satisfied and guarantee 
that the solutions of the series of approximations by AGMA 
converges to a point that satisfies the KKT conditions of fTOl 
i.e., a local maximum is attained. ■ 


Step 3: Stop if \\f3*{ti) — — 1)|| < ei, else set 

ti := ti + 1 and go to Step 1. 


C. Power Allocation Problem 


approximated monomial equalities and posynomial inequal¬ 
ities. Next, we show how we can transform the objective 
function into the monomial function to reach the GP based 
formulation for (fTol i. 

Proposition 2; Consider auxiliary variable xq > 0 and 
Si 3> 1. The user association problem (fTOl l at each iteration 
ti can be transformed into the following standard GP problem 


For a given set of UAFs derived from Algorithm A.I., the 
optimization problem of power allocation becomes 


max EE E E l^m,k,ng (j'')Rm,k,ng (P(^2)) 

^ mGAi g^Q UgGAfg kGK. 


( 20 ) 


subject to: 

C^2:EE E 

maM kCiK TLg^JSfg 

( 21 ) 


min xo(fi), (15) 

/3(tl).rCo(tl),Sm,ng(tl) 

subject to : C4.1-C4.4, 


yg&G, 

<^2:^ E <p^^ VmeA^, 

IT-gGAfg k^K. 


„ ra;o(U) 


-co(ii) 


^ l.Co(fl), 

l^m,k,ng {ti)R m,k,ng m) 


n 


mGA4,g^Q,ngGJ\fg ,kGfC 


< 1, 


(16) 


Cl.l : 


C-m^k^rig (f 1) 


n 

ng£j\fg,k£K *■ 


^m,k,ng (^1 


< 1 , 


C3.1; E E Pm,k,ngitl) < l, Vm e Ad,Vfc G /C, 

rag£A/g 


where t 2 is the index of iterations. Due to interference in 
the objective function of Rm,k,ng{P{t 2 )), (l20l) is a non- 
convex optimization problem. We again follow the approach of 
Section III. A to convert (l20l l into the GP optimization problem. 
First, we rewrite the objective of (l20l l as 

n Im.k,ng(Pit2)) (22) 

mGM,geG,ng£Arg,kCiK 


where 

7m,fe,ng (P(i2)) 
and 


Cr^ + Irr 


,(^ 2 ) 


At2)hr, 


Im^k,ngif2) 


where 

Pva^k^ngifl) — (12) 

I3m,k,ng{tl - l)Rm,k,ng{P{t)) 

SmeAt SrigeAAg Sfcg/C (^1 “ ^)Rm,k,ng{P{t)) 

yg&G. 

and Cm,k,ng{ti) and co(fi) are defined in (fTSI l and (fT^ at the 
top of the next page. 

Proof: See Appendix B. ■ 

Now, at each iteration, the optimization problem can be 
replaced by its GP approximation in (fTSl) . Iteratively, (fTSl l 
will be solved till the optimal value of UAF will be obtained. 
The overall iterative algorithm to solve (flSl l is summarized in 
Algorithm l.A where the iterative algorithm is repeated till 
the convergence condition holds. 

Proposition 3: With AGMA, Algorithm l.A converges to 


2m,fc,ng (fa) — 

E E E Pm',k,n'^ {t 2 )h m^k^n'g • 

Vn^e7Vg,ra^/ng 

Now from AGMA in Section III.B, 7m,fe,ng (P(f 2 )) can be 
approximated as 

^2 \ -«o(t2) 

7m,fc,ng (P(f2)) ~ “t” ^m^k^ng(j'2)) 


n 

mGM.,gGQ,ngGAfg ,kG}C 

where 


Pm,k,ngit2)hm,k,r. \ 

^m,k,ng (^ 2 ) 


^m,k,ng (^ 2 ) — 

Pm,k,ng (^2 l)^m,fc,ng 
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^m,k,ng (^1) — 

coih) = 


I3m,k,ng{ti - l)Rm,k,ngiP{t)) 


Xo{ti 1) + SgGS SnggTVg SfceA;l)^m,fc,ng (P(i)) 

Xoih - 1) 


Xo{ti 1) + SmGA< SgGS SrigGA^g SfeGK(^1 l)^in,fc,ng (P(0) 


(18) 

(19) 


Algorithm l.B: Algorithm to Derive Power Allocation 
Initialization: Set t 2 := 1 and obtain (3{t2 = 1) = (3{t) from 
Algorithm l.A. 

Repeat 

Step 1: Update Km,k,ngit 2 ) and Ko(< 2 ), 

Step 2: Find optimum power allocation according to (l23l l 
using CVX 111,11351, 

Step 3: Stop if ||P(f 2 )-P(f 2 -l)|| < £ 2 , else set t 2 := ^2 + 1 
and go to Step 1. 


and 

+ SmGAt.rigGA/'g.sGa P'rn,k,ng (^2 — ^)hm,k,ng 

Consequently, (l20l) is transformed into the following standard 
GP problem 

min 7 m,fc,ng(P(f 2 )) (23) 

^ mGA^,3GS,ngGA/’g ,fcG/C 

subject to: 

cr2: l[ 7m,fc,ng(P(t2))<2-^»\ 

m^A4,gGQ,ngGAfg,kGK, 

yg&G, 

<P^^ VmeM, 

g&G rigGAfg k^K 

The overall optimization problem is iteratively solved as 
described in Algorithm 1 .B until the power vector converges, 
i.e., ||P(f 2 ) — P(f 2 — 1)11 < £2 where 0 < £2 ^ 1- Note that 
Proposition III holds for the iterative Algorithm l.B. 


D. Computational Complexity 


When CVX is applied for the problems in Algorithms l.A 
and l.B, it solves GP problems with the interior point method. 
Given by ll35l . the number of required iterations is equal to 
, where c is the total number of constraints in (fTSl l. 
is the initial point for approximating the accuracy of interior 
point method, 0 < p <C 1 is the stopping criterion for interior 
point method, and ^ is used for updating the accuracy of 
interior point method ||35]| . In the problem described earlier, 
the number of constraints in (fTSl l is ci = G+MK+4MN+1 
for Algorithm l.A and C 2 = G + M for Algorithm l.B. 


Moreover, in Algorithms l.A and l.B, for each iteration, 
the number of computations required to convert the non- 
convex problems using AGMA into ( fTSl l and (|22]) is U = 
KM^N + 6KMN + MKGN and ^2 = GMKN + 2MKN, 
respectively. Therefore, the total number of computations for 


Algorithms l.A and l.B are 


n X 


Iog(ci/(t°gi)) 
iog(4i) ’ 


Computational complexity = < 


for Algorithm l.A, 

„• ,, log(c2/(t2g2)) 

*2 ^ Uite) ’ 

for Algorithm l.B. 


Based on this analysis, the computation complexity of 
Algorithm l.A is significantly higher than that of Algorithm 
l.B. Moreover, Algorithm l.A is more sensitive to K and N 
than Algorithm l.B. 


IV. Simulation results 
A. Simulation Parameters 

We consider a multi-cell VWN scenario with M = 4, 
K = A and G = 2. The radius of each BS is considered to be 
500 m where the users are randomly located (from a uniform 
distribution) within the whole area of interest unless otherwise 
stated. The channel gains are derived by the Rayleigh fading 
model in which hm,ng,k = Xm.ng.fcc(m^„g where /3 = 3 is 
the path loss exponent, dm,ng is the distance between the 
BS m and user Ug and Xm,n ,k is the exponential random 
variable with mean/parameter equal to 1. We set Ng = 4 
and P™™ = 40 dBm, unless otherwise stated. Also, for 
all of the simulations, we set 2i = 10^ and ei = 10“®, 
£2 = 10“®. In all of the following simulations, for each 
realization of network, when there exists no feasible solution 
for the system, i.e., C1-C4 cannot be hold simultaneously, the 
corresponding total rate is set to be zero. Also for all the 
following simulations, we have for all p € (/ and 

pmax _ p™ax all TO G Ad. To Compare the performance of 
our algorithm, we apply the traditional algorithm explained in 
the following subsection. 


B. Traditional Approach 

In practical scenario of wireless networks, e.g., 2G, the users 
are assigned to BSs based on the highest average received 
SINR measured during one measurement window. In this case, 
the resource allocation problem is transformed into 

E E E a ^'m.CngRm,k,ng{^) (24) 

m^M-g^Q Ug^Mg k^K 

subject to: Cl - C3, 

where /3' = [P'm,k,ng\ym,k,ng and shows the sub¬ 

carrier allocation of user Ug on sub-carrier k when it is allo¬ 
cated to the BS TO. To solve (l24l i. we can apply the algorithm 
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Fig. 3: versus (a) outage probability (Pr(outage)), and 

(b) total rate 


C. Evaluation of Algorithm 1 and Algorithm 2 


Fig. 2: Total rate versus (a) sub-carriers, K, and (b) 
maximum transmit power per BS, P™” 


similar to Algorithm 1 introduced in Section III, except that 
here, (fTsT i contains C1-C3, and C4 are removed. When the sub¬ 
carrier assignment is solved, the optimal power is derived from 
Algorithm l.B, for (l24l i. This iterative algorithm is terminated 
when the convergence conditions are met which is summarized 
in Algorithm 2. 


Algorithm 2 

Initialization: Set fa := 1, BS assignment where user Ug is 
assigned to BS m based on the average received SINR. 

Repeat 

Step 1: Compute /3 *(^ 3 ) from Algorithm l.B except that BS 
is assigned based on the signal strength. 

Step 2: For a fixed (3 *(^3), find the optimal power allocation 
P(f 3 ) from Algorithm l.B. 

Step 3: Stop if ||/3 *{h - 1)|| < £1 and ||P*(f 3 ) - 

P*(f 3 — 1)11 < £ 2 , where 0 < £2 <C 1 , otherwise := fa -f 1 
and go to Step 1. 


Primarily, we evaluate and compare the performances of 
Algorithm 1 and Algorithm 2 versus different number of sub¬ 
carriers and maximum transmit power in Figs. |2(a) and |2(b) 


respectively. The first point to notice from Figs. 2(a)| and 2(b)| 
is that Algorithm 1 considerably outperforms Algorithm 2 for 
different values of K and 

From Fig. |2(a)| it can be observed that the total rate is 
increased by increasing the total number of sub-carriers, i.e., 
K, due to the opportunistic nature of fading channels in 
wireless networks El. As also expected, with increasing 
the total achievable rate is also increased which is shown 
in Fig. |2(b)| Both figures indicate that by increasing the value 
of R^g'', the total rate decreases. It is because the feasibility 
region of resource allocation in dD is reduced by increasing 
R^g'', leading to less total achieved throughput ll35l . However, 
from Fig. |2(b)| increasing has considerable effect on 
the performance of Algorithm 2 compared to Algorithm 1. 
It can be interpreted as Algorithm 1 can efficiently control 
interference between different cells compared to Algorithm 2. 
Therefore, the chance of feasible power allocation for larger 
values of is increased by Algorithm 1. 

To study this point, we consider the outage probability of 
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Fig. 4: Illustration of network setup to investigate the 
coverage of multi-cell VWN 


Cl, expressed as 

Pr(outage) =Pr{ X! X! X! Pm^k^Ug Rm^k^rig (P? /3) 

mGM UgGAfg k£K 

Vpeg. 



(a) 



(b) 


Via Mont Carlo simulation, we derived Pr(outage) for the 
presented setup for both Algorithm 1 and Algorithm 2, de¬ 
picted in Fig. |3(a)| In this figure, we assume K = 8 and 
P™™ = 40 dBm for all m G M.. The results demonstrate 
that as the rate reservation per slice i?”'' increases, the outage 
probability for both Algorithm 1 and Algorithm 2 increases. 
However, Algorithm 2 has larger value of outage probability 
compared to the outage probability of Algorithm 1, meaning 
that the feasibility region of Algorithm 2 is smaller than that 
for Algorithm 1. On the other hand. Algorithm 1 can efficiently 
manage interference in the specihc region between different 
cells compared to Algorithm 2. It is mainly because Algorithm 
1 has more degree of freedom to choose the BS and allocate 
the sub-carriers among users of different slices. However, the 
BS assignment is predetermined in Algorithm 2. Therefore, 
the achieved rate derived by Algorithm 2 is less than that of 
Algorithm 1. Consequently, with increasing i?gthe decrease 
rate of Algorithm 2 is greater than that of Algorithm 1, since 
Algorithm 2 cannot manage the interference between different 
BSs. Hence, Algorithm 2 cannot satisfy the minimum required 
rate of slices leading to decrement of VWN efficiency. 

For the same setup, in Fig. |3(b)| the total throughput of 
Algorithms 1 and 2 are plotted for different values of P^''. Fig. 
|3(b)| clearly shows that Algorithm 1 yields higher throughput 
compared to Algorithm 2, while both of these algorithms be¬ 
come infeasible after = 4 bps/Hz. We should emphasize 
that in all simulation results, when the problem is infeasible, 
i.e., there is no power and sub-carrier vectors that can meet 
the constraint Cl for all g G G, the achieved total rate is set 
to zero. 


Fig. 5; Throughput distribution for (a) cell-edge, and (b) 
cell-center users 


D. Coverage Analysis 

In any cellular network, the coverage is one of the most 
important planning parameters which can be measured by 
SINR or achieved total rate of users at the cell boundaries. To 
study the performance of Algorithm 1 to increase the coverage, 
we consider the simulation setup similar to Fig. |4] where the 
majority of users are located in the cell-edge, consequently, 
these users experience high interference from other BSs. 
Therefore, the achieved rate of each user is decreased which 
can be considered as the worst-case scenario of coverage 
analysis. 

The cumulative distribution function (CDF) of the total 
throughput of cell-edge users and cell-center users are depicted 
in Figs. |5(a)| and |5(b)[ respectively for both Algorithms 1 and 
2. It can be seen that Algorithm 1 outperforms Algorithm 
2 for the cell-edge users where 50% of users in the cell- 
edge achieve a throughput of 2.5 bps/Hz by Algorithm 1, 
while their throughput is around 1.5 bps/Hz via Algorithm 
2. However, the performance of both algorithms are similar 
for the cell-center users. It is because via user-association in 
Algorithm 1, the interference among different cells can be 
controlled while Algorithm 2 cannot control the interference 
through the connectivity of users to different BS and it is per- 
determined by the received SINR of reference signal. In other 
words. Algorithm 1 can provide better coverage even for cell- 
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Fig. 6: Total rate for (a) uniform user distribution, and (b) 
non-uniform user distribution 

edge users for multi-cellular VWN which is desirable from 
implementation perspective. 

The performance is further investigated with respect to 
the number of users in the cell-edge in Figs. |6(a)| and |6(b)| 
Obviously, Algorithm 1 can consistently improve performance 
of cell-edge users and at the same time, maintain desirable 
throughput of each slice regardless of the user deployment 
density compared to the Algorithm 2. For instance, for N = 
18, for the uniform user distribution, the total rate is increased 
from 7 to 11 bps/Hz and from 24 to 32 bps/Hz for cell-edge 
and cell-center users from Algorithm 2 to Algorithm 1, i.e., 
57% and 33%, respectively. For the non-uniform case, when 
N = 32, the rate is incremented from 7 to 12 bps/Hz and 18 
to 27 bps/Hz for cell-edge and cell-center users, respectively, 
which is equal to 71% and 50% increment in total rate 
from Algorithm 2 to Algorithm 1, respectively. These results 
show the efficiency of applying Algorithm 1 in increasing the 
coverage over the whole network. 

E. Convergence Analysis of Algorithm 1 

In Fig. |7(a)| the number of iterations required for conver¬ 
gence for Algorithms 1 .A and 1 .B versus the total number of 
sub-caiTiers, K, is plotted for = 8 and = 2 bps/FIz. 
Similarly, in Fig. |7(b)| the number of iterations required for 
convergence versus the total number of users, N for K = A 


sub-carriers is plotted in the case of Algorithm LA and l.B. 
According to the computational complexity analysis in Section 
III.D, as N and K increase, the number of iterations required 
for convergence also increases. The computational complexity 
in the case of Algorithm 1. A is more than that of Algorithm 
1. B owing to the fact that the total number of constraints 
for Algorithm l.A is higher than Algorithm l.B which is also 
evident from the plots of iterations versus N and K in Fig. 

|7(i)|andFig.M 

V. Conclusion 

In this paper, we proposed the joint BS, sub-carrier and 
power allocation algorithm for multi-cell OFDMA based vir¬ 
tualized wireless networks (VWNs). In the proposed setup, we 
have considered a set of slices (service providers), each has a 
set of their own users and require a minimum reserved rate. We 
formulated the related optimization problem based on the new 
defined optimization variable, called user association factor 
(UAF) indicating the joint sub-carrier and BS assignment. 
To solve the proposed non-convex and NP-hard optimization 
problem, we followed an iterative approach where in each step, 
one set of optimization variables is derived. However, even for 
each step, the optimization problem is non-convex and NP- 
hard. To derive the efficient algorithm to solve them, we apply 
the framework of iterative successive convex approximation 
via complementary geometric programming (CGP) to trans¬ 
form the non-convex optimization problem into the convex 
one. Then, to efficiently derive the solution, we applied CVX 
to solve the optimization problem of each step in this called 
two-layer iterative approach. Our simulation results reveal that 
via the proposed approach, the throughput and coverage of 
VWN, specifically for the cell-edge users, are considerably 
improved compared to the traditional scenario where the BS 
is assigned based on the maximum value of SINK. 

Appendix 

A. Proof of Proposition 1 

From the definition of Xm,ng(ti) and yng(ti), C4 can be 
rewritten as 

Xm,ng{tl)[yng{tl) ~ Xm,ng{tl)] = 0 , (" 25 ) 

Vng e fig, Vp e Q, Vm e M, 

which is not a monomial function. (l25T l can be rewritten as 

Xni,ngiti)yng{ti) = and by adding 1 to both the 

left and right hand sides, we have Xm,ngiti)yngiti) + 1 = 
1 -I- xf^ ,^^{ti) for all ng € Afg, g € G, and m e M. We 
define Sm,n (ti) > 0 as an auxiliary variable to relax and 
convert (|25]) into the posynomial inequalities as follows ifTOl 

Xm,ng {tl)yng (ft) + 1 ^ Sm,ng (G) ^ 1 + 2^^ (^l)) (26) 

Vrig e Ng, Vg G G, Vto S M. 

The above inequalities can be written as 

Xm,ngifl)yngifl) + 1 ^ 

^m,ng (^1 ) 
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and 


^m^Ug (tl ) 


< 1 . 


1 +a;m.n,(il) 

Now, the above constraints can be approximated using AGMA 

approximation introduced in Section III. B. as 

C4.1: 4" ^rn.rtg (il)yng (^l) — 1) 


C4.2: [ 


1 




^7n,ng (G) 


r.^2 


2,ng (G) 


,(^l) 


, (*l) 


< 1 , 



_ O^m.ng (G ) _ 

where Xm,ng{tl) = ^ (t\-l) + l am,ng(tl) = 

n (*1-1) 

—^-^fjTTYypX- Now, C4 can be replaced by the following 
constraints 

C4.1: ■Sm,ng(^l) 4" (^l)yng (^l)'5m,ng (^l) — 

C4.2: (ti) 


(a) 


Xm,ng (fl) 

(i- 


^ ^m,n„ t''! 


—otf, 


,(*i) 


C4.3 

C44 


C^m,ng (fl)_ 

Xm,ng{tl) — ^ ^ / 
keK 


< 1, 




yng{tl) — ^ ^ Pm,k,ngiti), 

m^A4,k^K 

Note that via (l26l l. the positive condition for the con¬ 
straints of GP is met ||32l. However, the equality con¬ 
straints in C4.3 and C4.4 are not monomial since we 
have Xjn,ng{ti) - J2keK l^^,k,ngiti) = 0 and t/„g(fi) - 
T.m&M.keK l^rn,k,ng{ti) = 0, and, they have negative con¬ 
straints. To convert C4.3 and C4.4 to the monomial functions, 
we again apply AGMA approximation presented in Section 
III.A as 



C4.3: Xm,ngitl) 


keK 


C4.4: J] 

m^A4,kGK 


_ k^m,k,ng (^l) 


= 1 (27) 


k]m,k,ng (G) 


= 1 , 

(28) 

where Vm,k,ng{ti) and Pm,fe,„g(fi) are defined in ([T3]l and 
(Ell, respectively. 


(b) 

Fig. 7: Number of required iterations for lower-level iterative 
algorithms versus (a) sub-carriers, K, and (b) total number 
of users, N 


Clearly, we have negative elements on the objective function 
similar to our general formulation in Q. To meet the positive 
conditions of objective function in GP, we consider Si ^ 1 
and rewrite objective function as 

“1 ~ /^m.fc.rig (fl)77m,fc,7ig (P(f)) 

mGM g£0 ngGTVg fc£/C 

which is always positive. Then, considering a positive auxiliary 
variable Xq, and rewrite the objective function with this new 
auxiliary variables 

__ 

^0 + J2meM '^g&g SngGTVg 'l2keK Pm,k,ng (fl)7?m.fe,ng {P{t)) 

< 1, (29) 


B. Proof of Proposition 2 

To reach the GP based formula for (doll, we should have 
minimization over the objective function, i.e., 

min ^ ^ ^ —^m,k,ngRm,k,ng{P{t))- 

^ m^M gS:Q rig^Mg k^K 


Now, ( |29] ) can be rewritten as the product of monomial 
functions based on the AGMA from Section III. B as 


__ 

Xo 

^1 

.co(ii)_ 

\p 

Ti,k,ng 


* co(ii 


n 

T G A<, 9 e 5, Ti g G A/’g, fc G/C 

Cm,fe,ng(il) 


( 30 ) 


g (^l) 


< 1 , 
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where Cm,k,n (^i) and co(ti) are updated from (fTST i and ( fT9l l, 

respectively. Therefore, the total optimization problem can be 

transformed into (flSl l. 
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